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ABSTRACT 

For the study of deep inelastic scattering at HERA, the use of Jaquet-Blondel variables is 
advantageous in several respects. We calculate the complete leptonic O(a) QED corrections 
for the reaction ep —> eX in these variables. All but one phase space integrations are performed 
analytically. After exponentiation of soft photon corrections, an accuracy of at least 1% is 
matched. Numerical results are presented and compared to estimates based on the leading 
logarithmic approximation. 

Equation (17) is corrected in accordance with equation (7.44) of DESY 94-115. The 
misprint did not influence the numerics. 



1 Introduction 

At HERA, the deep-inelastic scattering of electrons off protons, 


e(ki) +p{ P i) -> e{k 2 ) + X(p 2 ), (1) 

allows the study both of the basic interactions of electrons and protons and of the proton 
structure in a completely new kinematic region [1], Reaction (1) has to be analyzed together 
with the corresponding photonic corrections from the radiative process, 

e(ki) +p(pi) -> e{k 2 ) + X(p 2 ) + n^{k). (2) 


Bremsstrahlung may substancially contribute to the observed net cross section. Its amount 
crucially depends both on the experimental cuts applied to the photonic phase space and on 
the kinematic variables and their allowed regions of variation. During a longer period, the use 
of leptonic variables dominated in the literature. At HERA, these have certain disadvantages: 
QED corrections may become extremely large, and in certain kinematic regions it is impossible 
to separate energetically the final electron from the bremsstrahlung photon. More advanta¬ 
geous are the hadronic variables on which the structure functions depend. Unfortunately, 
their experimental determination is also nontrivial or even impossible. A natural compromise 
has been proposed with the Jaquet-Blondel (JB) variables [3]. During the 1991 HERA Work¬ 
shop [2], features of the various sets of variables have been discussed in detail. Here we will 
concentrate on the case of JB variables. These variables are not influenced neither by the 
leptonic nor the longitudinal hadronic degrees of freedom: 
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with S = —{ki + Pi) 2 and Qh = P2 — Pi- Using the JB variables, one has the opportunity to 
search for a phase space parametrization with an early integration over the leptonic degrees of 
freedom. In another approach, we did not use this opportunity since we tried there to handle 
leptonic, mixed, and hadronic variables in a common formalism as much as possible [4, 5]. For 
details of the derivation of the formulae which will be presented here, we refer to that article. 

In this letter, we take advantage of the above-mentioned properties of the JB variables. We 
have performed all but one phase space integrations for the double-differential cross section. 
We explicitely scetch the rather compact complete O(a) leptonic corrections from reaction (2) 
to the double-differential cross section of reaction (1), including soft photon exponentiation. 
The leptonic corrections arise from the radiation of photons from the leptonic part of the Feyn¬ 
man diagrams. They are much larger than the other ones - quark- and lepton-quark interfer¬ 
ence bremsstrahlung. Further, leptonic bremsstrahlung may be treated model-independently , 
relying on some set of structure functions which need not necessarily be related to the quark- 
parton model. 

In Section 2, the kinematics is explained and the basic formulae are collected. Section 3 
contains numerical results and their discussion. 



2 Kinematics and phase space integrations 


2.1 The Born cross section 

The cross section of (1) is 
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where, of course, for the lowest-order cross section the different sets of kinematic variables 
agree; e.g. Q 2 B = Q 2 h , y JB = yh, Tjb = x h . Further, 

Sf = Q 2 , 

5 2 b = 2(1 -y)S\ 

5 3 b = 2(2 -y)Q 2 S. (5) 

The Ai are proportional to the structure functions 3 ( 37 , Ql) describing the electroweak 
interactions of leptons and nucleons; their arguments have to be defined by the hadronic 
kinematics: 
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where Qi, Vi and cii are the corresponding charge (Qi± = ±1), vector and axial-vector couplings 
of the lepton with the Z boson: vi = 1 — 41 Qi \ sin 2 6 W , ai = 1, and 9 W is the weak mixing 
angle. Further, 


x = x(Q 2 ) 


G E M 2 Q 2 

\/2 Sna Q 2 + M§ 


( 8 ) 


The F, G, H are structure functions parametrizing the general hadronic tensor for the |y| J , 7 Z, 
and \Z\ 2 contributions to the Born cross-section, respectively. 

For the radiative process, the sum at the right-hand side in (4) will become the integrand of 
a three-dimensional phase space integral. Further, the kinematic weights S t which accompany 
the structure functions become complicated expressions of the invariants. The specific problem 
which has to be solved for each set of variables separately is the subsequent analytic integration. 
In the next chapter, we comment on that in terms of Jaquet-Blondel (JB) variables and 
perform two integrations explicitcly. Then, we are left with a one-dimensional integral which 
has to be performed numerically. 



2.2 The phase space integration 

With JB variables, the following parametrization of the phase space may be derived 1 [5]: 

7 r /*Tmax 1 7 - — 777 ^ /*1 /*27T /* 

T = — J dQ] B dy iB J 2 dr——-— j ^dcos-d R ^ dip R j dY h . (9) 

Here, the 

rfr ^ = nS ( 5 ( 4 ) (P 2 -^Pi) ( 10 ) 

is the phase space element of the final hadron system which will be completely absorbed 
into the definition of the hadronic structure functions. Further, we use the invariant mass 
r = — (&2 + k ) 2 of the (e, 7 ) compound and the photonic angles d R , iy R in the rest system of 
this compound, defined by & 2 + k = 0. Further, m is the electron mass, and the upper limit 
of the r variation is 


Anax — (1 — 2 \Jb)(1 — UJb)S. 


In these variables, the double-differential cross section becomes: 
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( 12 ) 


It — Ui f 

Si(Q 2 jB,yjB,r) = ——J dcos‘d R d(p R S i (Q 2 j B ,yjB,T,cos'd R , <p R ). (13) 

The functions Si(Q 2 B , yj B , r, cos d R , tp R ) are calculated from the squared sum of the Feynman 
diagrams. Their explicit expressions may be found elsewhere [ 6 , 5]. In the derivation of (12) we 
took advantage of the fact that the hadronic variables, and thus the hadronic structure func¬ 
tions, are independent of the photonic angles. The hadronic transferred momentum squared 
Q\ and Q 2 B are related by the following relation containing the additional variable r: 

Ql = (P 2 ~ Pi) 2 = Qjb + V3B (t - m 2 ). (14) 

1 - vjb 

For Born kinematics, r approaches its lower limit and Q\ and Q 2 B agree. The integral over r 
in ( 12 ) has to be performed numerically as long as one cannot neglect the difference between 
Q 2 b and Q 2 - the structure functions are depending on Q 2 h and thus also on r. At the other 
hand, the integral over the photonic angles may be performed analytically. For details we 
refer to [7, 5] where the integration techniques are explained and most of the integrals used 
may be found. Again, we quote here only the final result which has been obtained by means 
of the program for analytical calculations SCH00NSCHIP [ 8 ]: 
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1 Whenever necessary, we retain the electron mass m although the ultra-relativistic limit is assumed. 
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2.3 The infra-red problem and soft photon exponentiation 

The integral (12) diverges at z 2 — 0 (k — 0). This may be seen from (15-17); the kinematic 
factors S{ contain the Born functions (5) together with the common, divergent factor (L r — 
2)/z 2 . The infra-red divergent part of (12) may be isolated: 
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From the above equations, it is easy to derive an explicit expression for the isolated infra-red 
divergent part d 2 a It contains an integral over r which may be performed explicitely using 
dimensional regularization [9, 5]. 

Finally, the resulting cross section is 
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As mentioned, from (22) one derives 

^rHQjbj 2/Jb) = y drF IR (QjB, Z/jb, t), (26) 

with F m defined in (23). The QED vertex correction <5 ver t compensates for the infra-red 
divergence V m (n) of the soft bremssstrahlung contribution, 

<5 V ert(QjB? 2/Jb) = — 2V m (n) (L JB — 1) — -Tjb^ + t^JB ^(l) — (27) 
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and Li 2 (l) is the Euler dilogarithm. The net photonic correction <5 VR is finite: 
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Finally, the effect of multiple soft photon emission may be taken into account by soft 
photon exponentiation: 
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3 Results and discussion 

The numerical results are obtained with an extended version of the FORTRAN program 
TERAD91 [5]. We use exactly the same specifications of kinematics, structure functions, weak 
Standard Model parameters for the Z amplitude, etc. as have been used at the 1992 HERA 
workshop [10]. For this reason, the running of the QED coupling has not been taken into 
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Table 1: Leptonic radiative corrections 5 of 0(a)in percent; first column this calculation, 
second column a LLA calculation [11], 


account in the matrix elements. We use first the 0(a) approximation. The table 1 contains 
the QED corrections 6, which are defined as follows: 
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The corrections are small for small Xjb and yj B . With rising yj B , they may become large. In 
magnitude, they fall in between the leptonic and the considerably smaller hadronic corrections. 

In the table, we also show numerical results from another calculation, which is restricted 
to the leading logarithmic approximation (LLA). These numbers have been taken from table 3 
in [10] and are based on formulae published in [11], The agreement is as good as one can 
expect for a LLA calculation. It is much better than has been observed for the case of 
leptonic variables, where certain kinematic singularities which are not related to the large 
logarithms from the particle masses, are much stronger pronounced. A naive conclusion from 
the comparison could be that experimentalists who apply radiative corrections to data should 
restrict to the use of the LLA formulae, thus optimizing the numerical computings. Let us 
remind here that this would not be too advantageous since both the exact O(a) and the 
approximated formulae are one-dimensional integrals. 




In figure 1, we show the same corrections, but now including soft photon exponentiation 
(SPE; broken curves). Its inclusion effects the cross sections only there, where the corrections 
are numerically large, i.e. near y — 1. 

To summarize , in this letter the complete 0(a) QED corrections, including soft photon 
exponentiation, to deep inelastic ep scattering are obtained in JB variables for the first time 
in form of a semi-analytical expression. The numerical improvement compared to a LLA cal¬ 
culation, which also contains a one-dimensional integration, is minor. 
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Figure 1: Leptonic radiative corrections 5 of 0(a ) without (solid curves) and with 
(broken curves) soft photon exponentiation in JB variables; parameter: xjb = 
0.0001,0.001,0.01,0.1,0.5. 




